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ABSTRACT 

Using a sample of 19,464 galaxies drawn from the DEEP2 Galaxy Redshift Survey, we study the 
relationship between galaxy color and environment at 0.4 < z < 1.35. We find that the fraction 
of galaxies on the red sequence depends strongly on local environment out to z > 1, being larger 
in regions of greater galaxy density. At all epochs probed, we also find a small population of red, 
morphologically early-type galaxies residing in regions of low measured overdensity. The observed 
correlations between the red fraction and local overdensity are highly significant, with the trend at 
z > 1 detected at a greater than 5a level. Over the entire redshift regime studied, we find that 
the color-density relation evolves continuously, with red galaxies more strongly favoring overdense 
regions at low z relative to their red-sequence counterparts at high redshift. At z > 1.3, the red 
fraction only weakly correlates with overdensity, implying that any color dependence to the clustering 
of ~ L* galaxies at that epoch must be small. Our findings add weight to existing evidence that 
the build-up of galaxies on the red sequence has occurred preferentially in overdense environments 
(i.e., galaxy groups) at z < 1.5. Furthermore, we identify the epoch (z ~ 2) at which typical ~ L* 
galaxies began quenching and moved onto the red sequence in significant number. The strength of the 
observed evolutionary trends at < z < 1.35 suggests that the correlations observed locally, such as 
the morphology-density and color-density relations, are the result of environment -driven mechanisms 
(i.e., "nurture") and do not appear to have been imprinted (by "nature") upon the galaxy population 
during their epoch of formation. 

Subject headings: galaxies:high-redshift, galaxies:evolution, galaxies: statistics, galaxies: funamental 
parameters, large-scale structure of universe 



1. INTRODUCTION 

The galaxy population both locally and out to z ~ 1 
is found to be effectively described as a combination 
of two distinct galaxy types: red, early-type galaxies 
lacking much star formation and b lue, late-type galax- 
ies with active star f ormation ( e .g.. IStrateva et al. 120011 
Baldrv et all l2004t iBell et all 12004 iMenanteau et all 
2006T ) . The spatial distribution of this bimodal galaxy 
population is frequently phrased today in terms of the so- 
called morpholog y-density or co l or-density relat i on. A s 
first maantified by jOemlerl (fl974l ). lDavis fc Gellerl (fl976T ) . 
and iDresslerl (|1980f ). the morphology-density relation 
holds that star-forming, disk-dominated galaxies tend 
to reside in regions of lower galaxy density relative to 
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those of red, elliptical galaxies. 

Many physical mechanisms that could be responsible 
for this correlation between galaxy morphology, star- 
for mation history, and environment have been proposed 
(see lCooper et alJl2006L for a review of probable mecha- 
nisms). Are the morphology-density and color-density 
relations a result of environment-driven evolution, or 
were these trends imprinted upon the galaxy population 
during their epoch of formation? Only through compre- 
hensive studies of galaxy properties and environments, 
both locally and at high redshift, will we be able to un- 
derstand the role of local density in determining the star- 
formation histories and morphologies of galaxies. 

While the close relationship between galaxy type 
and density was primarily uncovered via the study 
of nearby clusters, recent work using the 2-degree 
Field Galax y Redshift Survey (2dFGRS, iColless et all 
[200lL [200l and the Sloan Digital Sky Survey (SDSS, 
lYork et aLT l2000f ) has established that the connec- 
tions between local environment and galaxy proper- 
ties such as morphology, color, and luminosity extend 
over the f ull range of densities, from rich clusters to 
voids (e g., lKauffmann_gt^|2004j; IBalogh et all 12004; ' 
Blanton e t al.l 120051 : iCroton et al.l 120051 : iRoias et~al 



2005). Furthermore, high-resolution imaging and spec- 
troscopic data in increasingly more distant clusters 
(to z ~ 1) have shown that the trends observed lo- 
cally persist to hig her z, at least in the highest density 
environments ( e.g.. IBalogh et al.lll997h iTreu et al.|[2003: 
iPoggianti et al.ll2006h .~ 
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Using a large sample of ga laxies drawn from the 
DEEP2 Galaxy Redshift Survey. ICooperet alJ (|2006l ) ex- 
tended the understanding of the color-density relation at 
z ~ 1 across the full range of environments, from voids to 
rich groups, showing that the correlation between galaxy 
color and mean overdensity found locally is in place, at 
least in a global sense, when the universe was half its 
present age. While the role of environment appears to 
have been very critical at z ~ 1 and perhaps at earlier 
times, quantitative measures of the evolution of environ- 
mental influences on the galaxy population or of correla- 
tions such as the morphology-density and color-density 
relation are limited. 

Comparisons of local results with studies of high- 
redshift clusters have pointed towards significant evolu- 
tion in the relationship between galaxy pr operties and lo- 
cal environment from z ~ 1 to z ~ (e.g. , [Dressier et al.l 
[TM^lCouchlt~aIlfl99l: [Smith et al.ll2005l) . While these 
results indicate an environment-driven evolution in the 
galaxy population (i.e., pointing towards nurture versus 
nature as the origin of the galaxy bimodality) , such work 
has been limited to the vicinity of rich clusters, and thus 
we know little about evolution in the relationship be- 
tween galaxy type and density across the full scope of 
galaxy environments. Furthermore, clusters include only 
a relatively small fraction of the total galaxy population 
at any epoch by number (and an even smaller fraction 
by volume). Thus, the evolution of the color-density re- 
lation among the vast majority of the galaxy population 
remains unprobed. 

In large clusters, the physical mechanisms at work 
(e.g., galaxy harassment, ram-pressure stripping, and 
global tidal interactions) go beyond those acting in 
group-sized systems and the field. Results from the first 
comprehensive study of galaxy environment over a broad 
range of densities at high z indicate that such cluster- 
specific physical mechanisms cannot expl ain the global 
color -density relation as found at z ~ 1 (| Cooper et al.l 
2006). Accordingly, in looking for evolution in the rela- 
tionship between the bimodal nature of galaxy proper- 
ties and the local galaxy environment, we must turn our 
attention to the entire dynamic range of galaxy overden- 
sities at high redshift. 

In this vein, recent work employing a sample of low- 
and high-redshift fiel d galaxies from the V IM OS VLT 
Deep Survey fVVDS. ILeFevre et al.H2005bD has found 
a strong evolutionary trend in the color-density relation 
for galaxies spanning the redshift range 0.25 < z < 1.5, 
with the color-magnitude diagrams for galaxies at 0.9 < 
z < 1. 5 showing no signifi cant dependence on environ- 
ment (jCucciati et al J [2006) . Also working at high red- 
shift, a study of the blue fraction (that is, the fraction 
of galaxies that have bl ue color) in galaxy groups and in 
the field population bv lGerke et al.1 (|2006l ) finds instead 
that the field blue fraction significantly differs from that 
of the group population out to z ~ 1.3. 

In this paper, we use the large sample of high-z galax- 
ies obtained by the DEEP2 survey to conduct a detailed 
study of the color-density relation at 0.4 < z < 1.35. In 
§2, we discuss the data sample employed along with our 
measurements of galaxy environments and colors. Our 
main results regarding the relationship between color 
and environment are presented in §3. Finally, in §4 and 
§5, we discuss our findings alongside other recent results 



and summarize our conclusions. Throughout this pa- 
per, we assume a flat ACDM cosmology with fi TO = 0.3, 
Qa = 0.7, w = —1, and h = 1. 

2. THE DATA SAMPLE 
2.1. The DEEP2 Galaxy Redshift Survey 

The DEEP2 Galaxy Redshift Survey is a nearly com- 
pleted (> 95%) project designed to study the galaxy pop- 
ulation and large-scale structure at z ~ 1 (jDavis et al.l 
120031: iFaber et al J 120061) . To date, the survey has tar- 
geted ~ 50, 000 galaxies in the redshift range 0.2 < z < 
1.4, down to a limiting magnitude of Rab = 24.1. Cur- 
rently, the survey covers ~ 3 square degrees of sky over 
four widely separated fields. 

In this paper, we utilize a sample of 32,002 galaxies 
with acc urate redshifts (qu ality Q = 3 or Q = 4 as de- 
fined bv lFaber et~aTll2006l) in the range 0.4 < z < 1.35 
and drawn from all four of the DEEP2 survey fields. 
While at z > 0.75 the sample includes galaxies from each 
of th e four survey fields , spanning a total of > 10 point- 
ings (|Faber et al.l I2006T ). at 0.4 < z < 0.75 our sample 
selection is limited to galaxies in the Extended Groth 
Strip (EGS), wher e no color cut to p re-select for z > 0.7 
galaxies was used (jDavis et al.ll2006[ ). 

2.2. Measurements of Rest-frame Colors, Luminosities, 
and Environments 

Rest-frame (U — B) colors and absolute i?-band mag- 
nitudes, Mb, are calculat ed from CFHT B,R,I pho- 
tometry ([Coil et al.ll2"004bl) using t h e K-c orrection pro- 
cedure described in IWillmer et al.1 (|2006h . All magni- 
tudes di scussed within this paper are given in AB mag- 
nitudes |OkiXGun3[l983). For zero-point conversions 
between AB and Ve ga magnitudes, refer to Table 1 of 
IWillmer et all (p006h . 

For each galaxy in the data set, we compute the pro- 
jected 3 rd -nearest-neighbor surface density (S3) about 
the galaxy, where the surface density depends on the 
projected distance to the 3 rd -nearest-neighbor, -D p ,3, as 
S 3 = 3/(nDp 3 ). In computing S3, a velocity window of 
±1000 km/s is utilized to exclude foreground and back - 
ground galaxies. In the tests of ICooper et al.1 (2005), 
this environment estimator proved to be the most ro- 
bust indicator of local galaxy density for the DEEP2 
survey. To correct for the redshift dependence of the 
sampling rate of the DEEP2 survey, each surface density 
value is divided by the median S3 of galaxies at that 
redshift within a window of Az = 0.04.; correcting the 
measured surface densities in this manner converts the 
S3 values into measures of overdensity relative to the 
median density (given by the notation 1 + S3 here) and 
effect ively accounts for re dshift variations in the selection 
rate ([Cooper et al.ll2005l) . 

Finally, to minimize the effects of edges and holes in the 
survey geometry, we exclude all galaxies within 1 h" 1 co- 
moving Mpc of a survey boundary, reducing our sample 
from 32,002 to 19,464 galaxies. The redshift distribution 
of this subset is plotted in Figure Q] For complete de- 
tails regarding the computation o f the local environme nt 
measures, we direct the reader to lCooper et al.l ([2006). 

Within the DEEP2 sample, the bimodality of galaxy 
colors in rest-frame U — B color is clearly visible; out 
to z ~ 1.4, the galaxy color-magnitude diagram exhibits 
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a clear division into a relatively tight re d sequence and 
a mo re diffuse "blue cloud" of galaxies (W illmer et al.1 
2006) . For this study, we compute the fraction of galaxies 
on the red sequence using t he fol lowing color division, as 
defined by IWillmer et al.l ((2006) and illustrated as the 
dotted lines in Figure [2] and the dashed lines in the lower 
panels of Figure [3] 

U-B= -0.032(M B - 21.62) + G.454- 0.25 + 0.831. (1) 

The red fraction (/r) within a given redshift and envi- 
ronment range is given by the number of galaxies red- 
ward of this relation in U — B color divided by the total 
number of galaxies within the same bin of redshift and 
environment. 

As discussed by IWillmer et al.l ((2006), this division 
between the red sequence and blue cloud is derived 
from the work of Ivan Dokkum etall (|2000f ). in which 
the color-magnitude relation for early-type (i.e., solely 
morphologically-selected) galaxies is measured. The re- 
sults of this work and similar studies of clusters at 
z < 1 show that ellipticals form a red s equence in color- 
magnitude space with little scatter (e.g.. lEllis et al.llT997l ; 
IStanford et al.l 119981 ). Thus, dividing the sample ac- 
cording to the color-cut presented in Equation [1] effec- 
tively selects the red, early-type portion of the bimodal 
galaxy distribution. Details regarding the possible evo- 
lut ion of this colo r divis ion with redshift are highlighted 
by IWillmer et all (|2006l ) and examined in §3.1.1 of this 
work. 
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Fig. 1. — The observed redshift distribution for the 19,464 galax- 
ies in the four surveyed regions which are more than 1 h~ 1 comov- 
ing Mpc away from a survey edge and within 0.4 < z < 1.35. 
Galaxies at z > 0.75 are selected from all four DEEP2 fields, while 
at 0.4 < z < 0.75 galaxies are drawn from the Extended Groth 
Strip only; in that field, no color— cut is used to pre-select high- 
redshift (z > 0.7) sources. In the figure, a dashed line is employed 
to illustrate this transition in the sample at z = 0.75. The redshift 
histogram is plotted using a bin size of Az = 0.025. 



2.3. Sample Selection 

Because the DEEP2 survey extends over a broad red- 
shift range, selecting galaxies according to a fixed ap- 
parent magnitude limit introduces differences in selec- 
tion with redshift that depend in a nonnegligible way 
on galaxy color and luminosity. To study how these se- 
lection effects potentially influence our results, we iden- 
tify several subsamples within the full catalog (Sample 
A) of 19,464 DEEP2 galaxies spread across the redshift 



M cut (zM-B) 



range 0.4 < z < 1.35 (cf. Fig. 1). The simplest selection 
method is to produce a volume-limited subsample ac- 
cording to a strict cut in absolute magnitude. We create 
such a sample (Sample B) by restricting to 0.4 < z < 1.1 
and requiring Mb < —21, the absolute magnitude to 
which DEEP 2 is complete along both the red sequence 
and the blue cloud at z = 1.1 (cf. top of Figure [5]). 

Producing a volume-limited sample with a fixed 
absolute-magnitude limit at all colors severely restricts 
either the redshift range probed or the number of galax- 
ies selected at each redshift. Here, we adopt a limiting 
redshift of z = 1.1 in order to maintain enough galaxies 
(2,784 in Sample B) with which to accurately compute 
the red fraction over the full range of overdensities. Given 
the large redshift range probed in our analysis, however, 
the number of galaxies with Mb < —21 at low redshift 
(0.4 < z < 0.75) is fairly small, and as such results at 
low z using Sampl e B are quite noisy. 

As discussed by Ge rke et al.l (|2006l ). in studying the 
evolution of galaxy properties it is also possible to pro- 
duce volume-limited catalogs with a color-dependent, 
absolute-magnitude cut by defining a region of rest- 
frame color-magnitude space that is uniformly sampled 
by the survey at all redshifts of interest. For the DEEP2 
survey, such a selection cut is illustrat ed in the bottom o f 
Figure[2]of this paper and Figure 2 of lGerke et alj (2006) 
and given by 

Q(Z - Zlhn) + 

mm {[a(U -B)+b], [c(U - B) + d]} , 

(2) 

where zy lm is the limiting redshift beyond which the se- 
lected sample becomes incomplete, a, b, c and d are 
constants that are determined by the limit of the color- 
magnitude distribution of the sample with redshift z < 
zii m , and Q is a constant that allows for linear red- 
shift evolution of the typical galaxy absolut e magnitude 
Mr. For the parameter Q, we adopt the iFaber et~aLl 
(2005) value of Q — —1.37, determined from a study of 
the i?-band galaxy luminosity fu nction in the COM BO- 
17 (|Wolf et al.l [200lh . DEEP1 dVotrt et all 12005 and 
DEEP2 surveys. 

By including this linear M B evolution in our selection 
cut, we are selecting a similar population of galaxies with 
respect to M B at all redshifts. Adopting this approach 
with a limiting redshift of zu m — 1.3, we define a sample 
of 11,192 galaxies (Sample C) over the redshift range 
0.4 < z < 1.3 that is volume-limited relative to M B and 
selected according to this color-dependent cut in Mb- 
The values of the constants a, b, c, and d which define the 
color-dependent selection are —2.06, —18.9, —2.63, and 
— 18.5, respectively. For compl ete details of the se lection 
method, we refer the reader to I Gerke et al.l (|2006T ) . 

As a final sample, we again select galaxies assuming 
linear evolution in M B , as for sample C. However, in de- 
termining the limiting absolute magnitude as a function 
of redshift, we do not apply a color-dependent selection 
as utilized for Sample C. Instead, we restrict galaxies of 
all color to the selection limit, 

M B {z) < M lim - Q(z lim - z) 

or (3) 
M B (z) - M*(z) < M lim - Af*(z lim ), 

where the limiting redshift is chosen to be zn m = 1.1 and 
M\i m = — 21 defines the absolute magnitude at z = 1.1 to 
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Fig. 2. — (Top) The rest— frame color— magnitude diagram for 
galaxies in the redshift range 1.075 < z < 1.125, selected from 
the full galaxy catalog (Sample A). The dotted vertical line defines 
the color— independent completeness limit of the DEEP2 survey 
at this redshift — the absolute-magnitude limit used for Samples 
B &; D. To this limit, the samples are complete for galaxies of 
all color at z < 1.1 (Bottom) The rest-frame color-magnitude 
diagram for galaxies in the redshift range 1.275 < z < 1.325, drawn 
from the full galaxy catalog (Sample A). The dotted line defines the 
completeness limit of the DEEP2 survey as a function of rest— frame 
color at redshift z = 1.3, as given by Equation [2] Due to the R— 
band magnitude limit of the survey, the fraction of red galaxies in 
the overall DEEP2 sample decreases significantly at high redshift, 
but by staying within this limit we can construct samples that are 
free of this effect. The dashed horizontal line shown in both plots 
illustrates the division between the red sequence and the blue cloud 
used t hroughout this paper, following the relation of Will mer et al.l 
(2006) (cf. EquationQ}. 

which DEEP2 is complete along both the red sequence 
and the blue cloud (cf. the top panel of Fig. [5]). Due to 
the more severe, color-independent cut in absolute mag- 
nitude, Sample D totals only 2,150 galaxies, spanning the 
redshift range 0.4 < z < 1.1. In contrast, while Sample 
B is similarly selected based just on absolute magnitude 
(that is, independent of galaxy color), the magnitude cut 
in Sample D evolves with redshift, staying fixed relative 
to Mg. Thus, Sample B selects galaxies down to the 
same absolute magnitude at all redshifts, and Sample 
D samples galaxies to the same depth in the luminosity 
function at each z. A brief summary of the four galaxy 
samples utilized in this paper is provided in Table [TJ 

2.4. Mock DEEP2 Survey Catalogs 

In order to test for possible systematic effects we em- 
pl oy a set of 12 m ock galaxy catalogs based on those 
of lYan et al.l (|2004h . These catalogs are derived from 
iV-body simulations by populating dark matter halos 
with galaxies according to a halo occupation distri- 
butio n (HOD) function (jPeacock fc Smith! 120001: iSeliakl 
2000), which describes the probability distribution of 
the number of galaxies in a halo as a function of the 



host halo mass. The luminosities of galaxies are then 
assigned according to the condition al luminosity func - 
tion (CLF) formalism introduced by lYang et alj (2003), 
which allows the galaxy luminosity function to be mass- 
dependent as well. Parameters for the HOD and the 
CLF are chosen to match th e 2dFGRS luminosity func- 
tion (iMadgwick et al.l I2002D an d two-point correlation 
function (|Madwgick et a l. 2003). By assuming that the 
manner in which dark matter halos are popu lated with 
galax ies does not evolve from z ~ 1 to z ~ (Ya n et ahl 
|2003|) save via an overall evolution in (correspond- 
ing to Q — —1 here), mock catalogs can then be built 
using dark-matter-only simulation outputs at varying 
redshifts. The re sulting simulated galaxy catalogs from 
lYan et al.l (pOOl are in excellent agreement with the 
lower redshift (0.7 < z < 0.9) DEEP2 correlation func- 
tion (ICoil et alJl2004bD an d the COMBO-17 luminosity 
function (|Wolf et alJl2003h . 

It is critical for this paper that we characterize any 
systematic effects which alter the relationship between 
galaxy color and observed overdensity as a function of 
red shift. The r efore, we employ a modified version of 
the lYan et all 



mock catalogs which incorporate 
galaxy colors; the se catalogs are described in detail by 
iGerke et al.l (|2006h . 

To construct these new catalogs, we first measure the 
local overdensity of each object in the lYan et all (|2004l ) 
mock catalogs using the full, volume-limited catalog. 
Each galaxy is assigned a rest-frame U — B color drawn 
from the set of real DEEP 2 galaxies at 0.8 < z < 1.0 
located in a corresponding bin of local galaxy overden- 
sity and absolute £?-band m agnitude, emp l oying the en- 
vironment measurements of ICooper et al.l ((2006). This 
produces mock catalogs which reproduce the observed 
relationships between galaxy color, luminosity, and envi- 
ronment in the DEEP2 data at 0.8 < z < 1.0, but do not 
allow those relationships to evolve with redshift unless 
the "true" distribution of galaxy overdensities evolves. 
After galaxy colors are assigned, apparent R magnitudes 
are determin ed using the K-corre ction methods applied 
for DEEP2 (jWillmer et al.l 120061 ). so that red galaxies 
will be lost at the same luminosity at a given redshift as 
in the data. We then apply the standard DEE P2 target- 
selection and slitmask- making procedures ([Davis et al.1 
120031: iFaber et al.l l2006) to these mock catalogs, so that 
we may directly determine the effects of DEEP2 target- 
selection algorithms on observed trends. 

These mock catalogs should not be a perfect repre- 
sentation of reality, as uncertainties in the observed en- 
vironments will cause some objects to be assigned the 
color of a galaxy that is actually in a higher- or lower- 
density environment than was measured; however, this 
effect is small (i.e., environment measurement errors at 
z = 0.8 — 1 are significantly smaller than the bin sizes). 
Nevertheless, they provide a robust test of whether we 
may falsely observe an evolution in the relationship be- 
tween galaxy color and environment when in fact none 
exists. They also allow us to directly explore what ef- 
fects the increased errors in environment measurements 
at higher redshifts (where samples grow dilute) may have 
on our results. We discuss these tests in S3. 2. 



3. RESULTS 
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TABLE 1 

Descriptions of Galaxy Samples 



Sample 


-^galaxies 


z range 


Brief Description 


Sample A 
Sample B 

Sample C 
Sample D 


19,464 
2,784 

11,192 
2,150 


0.4 < z < 1.35 
0.4 < z < 1.1 

0.4 < z < 1.3 
0.4 < z < 1.1 


all galaxies after boundary cut 

color-independent, volume-limited (Mg < —21) cut 
color-dependent limit, with limit held constant rela- 
tive to Mg (z) 

color-independent limit, with limit held constant rel- 
ative to M B (z) and set as Mb = —21 at z = 1.1 



Note. — We list each galaxy sample employed in the analysis, detailing the selection cut used to define the sample as well as the number 
of galaxies (JVg a i ax i cs ) included and the redshift range covered by each sample. 



3.1. The Color-Density Relation at 0.4 < z < 1.35 

Figure [3] shows the relationship between the red frac- 
tion (/ R ) and local galaxy overdensity for the full data 
set (Sample A) in four distinct redshift bins spanning 
0.4 < z < 1.35. We find that the red fraction exhibits a 
strong dependence on local environment in each redshift 
bin, such that /r is higher in regions of greater overden- 
sity. Fitting a linear relation to the trends in Figure 03 
we find that the slope of this color-environment corre- 
lation evolves with z, with the relative fraction of red 
galaxies in dense environments decreasing with lookback 
time (cf. Table [2]) . Over the full redshift range probed 
by the DEEP2 sample, the measured slope decreases 
from d/ R /dlog 10 (l + S 3 ) = 0.111 for 0.4 < z < 0.75 to 
d/ R /dlog 10 (l + 6 3 ) = 0.046 for 1.0 < z < 1.35, with the 
difference being significant at a greater than 3-ct level. 

The vertical error bars in Figure [3] indicate the Pois- 
sonian uncertainty in each point. Error estimates based 
on bootstrap and jackknife resampling amongst the 10 
DEEP2 pointings used yielded comparable uncertainties 
(within 10-20%), suggesting that sample (or "cosmic") 
variance is not the dominant source of uncertainty; while 
cosmic variance will influence the distribution of envi- 
ronments in a given redshift bin, it should not strongly 
affect the relationship between galaxy color and density 
in a given environment bin. Refer to §4.1 for further in- 
vestigation of the role of cosmic variance in this analysis. 

All of the samples described in §2.3 yield similar results 
for the evolution of the color-density relation (cf. Fig. [4] 
and Table [2]) . While overall /r, values by overdensity dif- 
fer when either a color-dependent or a color-independent 
Mb limit is applied to the full sample, in every case the 
color-environment relation weakens, but still is present, 
at higher z. Due to their small sample sizes, the effect is 
not statistically sigificant for Samples B and D, however. 

To further illustrate the evolution in the color-density 
relation with z, we present color-magnitude diagrams in 
Figure [5] divided by local environment and redshift. In 
this figure, the relationship between red fraction (/r) 
and galaxy overdensity is made directly discernable to 
the eye; in high-density regions, there is an increased 
proportion of sources on the red sequence relative to that 
observed in low-density environments. 

While the results shown in Figure [3] show a clear sig- 
nature of evolution in the color-density relation over the 
redshift range 0.4 < z < 1.35, relatively large redshift 
bins are used, thereby coarsely sampling the redshift do- 
main. To study the evolution of the /R-overdensity re- 
lationship with z in more detail, we divide Sample C 
into thirds according to overdensity and compute the 



red fraction as a function of redshift in a sliding bin 
of width Az = 0.1 for the galaxies in the high-density 
(log 10 (1 + S 3 ) > 0.25) and low-density (log 10 (1 + 5 3 ) < 
—0.2) extremes of the overdensity distribution, with re- 
sults shown in Figured Here, the high-density and low- 
density thirds are defined with respect to the overdensity 
distribution in each z bin; however, dividing the galaxy 
sample into high- and low-density regimes according to 
the overdensity distribution in the entire redshift range 
(0.75 < z < 1.3) probed yields consistent results. As 
detailed in §2.3, sample C is defined so as to be com- 
plete with respect to M B across the full redshift range, 
0.4 < z < 1.3, establishing a sample uniformly selected 
at all redshifts studied. 

As illustrated in Fig. [6j the color-density relation 
shows a continuous evolution with redshift from z ~ 0.75 
to z ~ 1.3, such that at redshifts approaching z ~ 1.3 
the red fraction in low- and high-density regions are sta- 
tistically consistent with each other. Extrapolating lin- 
ear regression fits to the /r(z) relations in high-density 
and low-density environments, we find convergence at a 
redshift of z = 1.32. We investigate this convergenece 
in more detail in the following subsections (§3.1.2 and 
§3.1.3). 

In our analyses, we have attempted to carefully account 
for selection effects due to the design parameters of the 
DEEP2 survey. For all galaxy samples but Sample A, se- 
lection effects related to the survey's magnitude limit and 
i?-band selection, which bias the sample against faint 
and/or red galaxies at higher redshifts, shoud be mini- 
mal. The agreement amongst the highly disparate sam- 
ples indicates that the results presented above are robust 
to such effects. 

The se results are similar to the findings of iGerke et al.l 
(2006), which investigates the evolution of the blue frac- 
tion among group and field galaxies at 0.7 < z < 1.3 
in DEEP 2. That work employed a very different tech- 
nique for measuring galaxy environment and studied a 
somewhat different set of galaxy subsamples. The mea- 
su rements of the evolu tion in the blue fraction presented 
bv lGerke et al.l (|2006h . however, are significantly noisier 
than those presented in this study and could be suscep- 
tible to biases if the group-finding algorithm used per- 
forms differently at different z. We note that while that 
work uses a different methodology for dividing galaxies 
into environment bins (identifying members of galaxy 
groups rather than with a continuou s meas ure of local 
galaxy density), both the lGerke et al.l (|2006[ ) results and 
those presented here are derived from subsamples of the 
DEEP2 survey; hence any systematics affecting one may 
also affect the other. 
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TABLE 2 

Linear Fits to the /r-overdensity Relation 
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Note. — We list the coefficients and 1— it uncertainties for the parameters of the linear-regression fits to the red fraction versus median 
overdensity relation given by /r = ao ■ log 10 (l + £3) + a\ (cf. Fig.[3]and Fig. 2} in each 2 bin employed for all four galaxy samples used. We 
also give the number of red— sequence members (AT re( j) along with the total number of galaxies (JVtot) in each redshift bin for that sample. 
For details regarding the various galaxy samples, refer to §2.3 and Table [T] 



3.1.1. Effect of Possible Evolution in the Color of the 
Bimodality 

While the position of the bimodality in U — B versus 
Mb color-magnitude space shows no significant depen- 
dence on redshift within the DEEP2 data, when studying 
the red fraction as a function of redshift, we must con- 
sider the possibility that there is evolution in the color of 
the red-sequence population and thus in the position of 
the color bimodality. If the location of the color bimodal- 
ity evolves with reshift, then the measured evolution in 
the red fraction in high-density environments could sim- 
ply result from galaxies moving (e.g., via passive evolu- 
tion) across our non-evolving color division (cf. Equation 
[T]). However, any affect related to possible evolution in 
the location of the U — B bimodality would be minimized 
by the differential nature of our measurements. That is, 
the location of the color bimodality does not appear to 
depend on environment, and thus any evolution in the 
red fraction for the galaxy sample in high-density re- 
gions should be mirrored in the /r(z) for the sample of 
galaxies in intermediate- and low-density environments. 

Furthermore, given that passive evolution will cause 
galaxies to move redward as they evolve and that 
the U — B color division in Equ ation [1] is defined 
by a galaxy sample at z ~ 0.8 (jWillmer et al.l 120061 : 
Ivan Dokkum fe Franxl I2001D , a non-evolving color di- 
vision will not increase the measured evolution in the 
color-density relation due to an interloper population of 
blue galaxies on the red sequence at z > 1. Instead, 
a non-evolving cut would create a more stringent selec- 
tion criterion at higher redshifts, limiting the number of 
galaxies residing in the trough between the red sequence 
and blue cloud (also called the "green valley") that are 
counted as members of the red sequence, thereby actually 
producing weaker evolution in the color-density relation. 

In the spirit of thoroughness, however, we test the ef- 
fect of using an evolving color division on our results. Al- 



lowing for passive evolution in the location of the color bi- 
modality of ~ 0.15 magnitudes (in U—B) p er unit z (e.g., 
Ivan Dokkum fc Franxj|200lHBlantonll2006h . such that the 
division between red and blue galaxies moves blueward 
with lookback time, we find no significant change in our 
results. That is, employing a redshift-dependent division 
between the blue cloud and red sequence, the relative 
red fraction in high- and low-density environments still 
shows strong evolution, consistent with no color-density 
relation at z ~ 1.35. 

3.1.2. Effect of Redshift Incompleteness 

In analyzing the evolution of the red fraction (/r ) as a 
function of redshift, we must consider the possible impact 
of redshift-dependent selection effects in the DEEP2 sur- 
vey. In particular, due to the finite amount of slit-length 
real estate available on DEIMOS slitmasks and the finite 
amount of observing time dedicated to the project, the 
DEEP2 survey only spectroscopically observes ~ 60% 
of the galaxies that meet the survey's target-selection 
criteria. Furthermore, about 30% of those galaxies tar- 
geted for spectroscopy fail to yield a high-quality (quality 
Q = 3,4) redshift. 

Initial follow-up observations of sources for which 
DEEP2 fails to measure a redshift indicate that roughly 
half of all failures (i.e., ~ 15% of DEEP2 targets) are 
at redshifts beyond the range probed by DEEP2 (i.e., 
z > 1.4; C. Steidel, private communication). Addition- 
ally, redshift failures may result from poor observing 
conditions, data reductions errors, or instrumental ef- 
fects. Clearly important for this work is the fact that the 
DEEP2 redshift failure rate is correlated with observed 
galaxy color and magnitude, thus possibly introducing 
bogus trends in the measured Jr(z) relation. 

To understand the degree to which our results are im- 
pacted by redshift incom pleteness, we employ t he weight- 
ing scheme presented by IWillmer et all {2006). The de- 
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Fig. 3. — Top Row: Using the full galaxy sample (Sample A), we plot the red fraction as a function of overdensity, log 10 (l + 53), in 
four distinct redshift ranges (redshift limits indicated at the top of each column). The circular points give the red fraction as function of 
the median overdensity, computed in four distinct bins of log 10 (l + £3). The horizontal error bars run from the twenty— fifth percentile 
to the seventy-fifth percentile of the overdensity distribution in each bin. The vertical error bars give the 1—a uncertainty on the red 
fraction within each overdensity bin, given by Poisson statistics. The grey shaded region in each panel shows the l—tj range of the red 
fraction as a function of median overdensity in a sliding bin with width given by the horizontal line segment in the upper left of each plot 
(A log 10 (1 + 53) = 0.2). The number in the upper left corner denotes the total number of red galaxies in the particular redshift interval. 
The dashed line in each panel shows a linear-regression fit to the four data points, with coefficients given in Table [2] Bottom Row: We 
plot the rest-frame color— magnitude relation for all objects in each r edshif t bin. The division between the red sequence and the blue cloud 
is given by the dashed line, following the relation of Willmcr ct al. (2006) (cf. Equation [TJ • The total number of galaxies (both blue and 
red) in each z bin is enumerated in the bottom right corner and given in Table [2] 
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rived weights account for variations in both redshift in- 
completeness and targeting rate as a function of ap- 
parent R magnitude and apparent R — I and B — R 
colors, assuming that the redshift distribution for red 
(R — I > 1.03) failures is identical to that of the ob- 
served (quality Q = 3,4) red galaxy sample and that 
the blue (R — I < 1.03) failures sit beyond the redshift 
range of the survey (defined as the "optimal" model in 
IWillmer et all 12000) . If we use these weights to com- 
pute fn(z) for Sample C (as in Fig. 0, all changes in 
the measured trends with redshift are well within the 
measurement uncertainties. Thus, biases related to red- 
shift incompleteness should not introduce spurious re- 
sults in our analysis. With and without weighting for 
incompleteness, we find a convergence of the red frac- 
tion in high- and low-density environments at z > 1.3; 
extrapolating linear regression fits to the fn(z) relations 
in high-density and low-density environments, we find 
convergence at a redshift of z = 1.33 when weighting for 
incompleteness versus a convergence at z = 1.32 without 
weighting. 

3.1.3. Tests with Mock Galaxy Catalogs 

The mock galaxy catalogs described in §2.4 allow us 
to further test the degree to which our observations of 
the color-density relation may be subject to systematic 
effects within DEEP2 (e.g., due to increased errors in en- 
vironment measures as the sampling density drops with z 
or due to the expected changes in the magnitude of pecu- 
liar velocities with z compared to the 1000 km/s window 
used). These catalogs should exhibit no evolution in the 
color-density relation when measured with perfect infor- 
mation, by construction. Therefore, when we apply our 
measurement techniques to subsets of the mock catalogs 
which emulate the observed samples, any apparent evo- 
lution in the color-density relation will generally indicate 
the presence of observational biases. 

There is a notable exception to this. It is possible to 
select a subsample from the volume-limited mock galaxy 
catalog (i.e., the catalog before applying DEEP2 target- 
selection criteria or slitmask-making algorithms) us- 
ing a color-dependent, absolute-magnitude limit which 
evolves with z, as for our Sample C above; however, for 
the mock catalogs, Q = — 1 is the appropriate evolution, 
as that value was used in their construction, rather than 
Q = —1.37 as used for Sample C. We find that if we use 
the full volume-limited mock catalog to determine local 
overdensities in real-space, there is in fact evolution in 
the color-density relation for such a sample (cf. Figure 

Et). 

This evolution does not reflect an observational se- 
lection effect, however. Instead, it is a result of ap- 
plying an absolute-magnitude cut which changes with 
z to mock catalogs which were constructed assuming 
rest-frame color depends on Mb in a way which is in- 
dependent of z (cf. §2.4). Because M B was brighter at 
higher z, the galaxies included by DEEP2 tend to be 
brighter at higher Z] and as may be seen in Fig. [2] and 
Fig. [51 brighter galaxies are more likely to be red. The 
color-dependent, absolute-magnitude limit exaggerates 
this effect. If instead we apply a color- and redshift- 
independent absolute-magnitude limit to the volume 
limited mock sample, the mock catalogs show consider- 
ably weaker color-density evolution, as we would expect 



(cf. Fig. [7p) . For panels (a) and (b) of Figure [3 we de- 
termine the "true" local galaxy environment according 
to the 7 th -nearest-neighbor surface density (converted 
into an overdensity, (1 + £7)) using the real-space galaxy 
positions and the full volume-limited galaxy population. 

Although the mock catalogs exhibit a modest evolu- 
tion in the color-density relation for samples with an 
absolute-magnitude limit that shifts with M* (like our 
Sample C), they still may be used to determine the net 
effects of DEEP2 sample selection on the observed color- 
density relation. We do this by comparing fn(z) in two 
different mock catalog samples. For the first, we do not 
apply DEEP2 target-selection and slitmask-making cri- 
teria, but only the color-dependent absolute-magnitude 
limit given by Equation [2] with Q — —1, and we use 
the best possible measure of environment in the cata- 
logs (based upon the real-space overdensity of galax- 
ies in the full, volu me-limited catalog, as employed by 
iCooper et al.l 12005) to divide the sample into environ- 
ment bins; see Fig. [7^. For the second, we apply the 
DEEP2 target-selection and slitmask design algorithm 
to these mock catalogs, and use the overdensity measured 
from this "observed" sample to partition the sample; re- 
sults for this sample are shown in Fig. [7fc. 

We find that in fact the DEEP2-like sample exhibits 
both a smaller gap between /r(z) for the extreme en- 
vironment bins and a stronger apparent evolution over 
0.75 < z < 1.3. This is consistent with a scenario where 
noise in environment measures causes objects to cross 
the boundary between the intermediate environment bin 
and the extremes, thereby artificially increasing the red 
fraction for the low-density sample while decreasing /r 
for the high-density sample. This effect should become 
worse at higher 2, as the sampling density decreases 
there, leading to greater environment errors. Therefore, 
though we observe convergence of between the red frac- 
tion in low- and high-density environments at z = 1.32, 
in actuality convergence should not occur until a some- 
what higher redshift. 

However, this test shows that observational selection 
effects cannot explain the convergence in /r that we ob- 
serve in high-density and low-density regions at z ~ 1.3. 
Extrapolating linear-regression fits to the /r(z) relations 
in high-density and low-density environments, we find 
convergence at z — 2.50 if we measure environment with 
the full sample (Fig. [7^,), or z = 2.29 after we apply 
DEEP2 target selection procedures (Fig.[7h). The mock 
catalogs provide clear evidence that the null hypothe- 
sis of no evolution in the color-magnitude relation over 
0.75 < z < 1.3 cannot hold. 

3.2. Red Galaxies in Low-Density Environments 

The relations presented in Figure [3] also show evi- 
dence of a trend with z in the normalization of the /r- 
overdensity relation, such that the total fraction of red- 
sequence galaxies averaged over all densities decreases 
with z. This effect is dominated by two systematic trends 
in the data: (1) from z ~ to z ~ 1, the number density 
of red gala xies in the universe decreases by a factor of 
~2-3 (e.g.. iBell et all[200l iFaber et al.ll2005h . and (2) 
at redshifts beyond z ~ 1.1, the red galaxy fraction de- 
creases precipitously in the DEEP 2 sample, due to th e 
survey's i?-band magnitude limit ([Cooper et alj [20061 . 
This latter effect is quite apparent in the final redshift 
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Fig. 4. — We plot the red fraction as a function of overdensity, log 10 (l + 53), in distinct redshift ranges for the remaining three galaxy 
samples identified in §2.3 and summarized in Table[T] — Samples B (Top Row), C (Middle Row), and D (Bottom Row). The circular points 
give the red fraction as function of the median overdensity, computed in distinct bins of log 10 (l + S3). The horizontal error bars run from 
the twenty-fifth percentile to the seventy— fifth percentile of the overdensity distribution in each bin. The vertical error bars give the 1— a 
uncertainty on the red fraction within each overdensity bin, given by Poisson statistics. The grey shaded region in each panel shows the 1— a 
range of the red fraction as a function of median overdensity in a sliding bin with width given by the horizontal line segment in the lower 
right of each plot (Alog 10 (1 + £3) = 0.2). The number in the upper left corner denotes the total number of red galaxies in the particular 
redshift interval. The dashed line in each panel shows a linear— regression fit to the data points, with coefficients of the fit given in Table 
[2] For the each sample, the baseline value of /j^ varies significantly as the different sample selection criteria include different portions of 
color-magnitude space. 
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Fig. 5. — We plot the color-magnitude relation for DEEP2 galax- 
ies in low-density and high— density environments at low and high 
redshift. Here we use "low" and "high" density to refer to the 
extreme thirds of the overall overdensity distribution within Sam- 
ple C. The dashed line in each plot denotes the division between 
the red sequence and blue cloud as given in Equation [T] At the 
bottom right corner of each plot, the red fraction (/r) for that 
color-magnitude diagram is given. For each row, the two panels 
contain equal numbers of objects. The greater prevalence of red 
objects in high-density regions, particularly at lower redshifts, is 
readily apparent to the eye. 

bin in Figure [31 where the normalization of the observed 
red fraction-overdensity relation drops significantly; the 
other volume-limited samples studied here should not be 
affected by this, however (cf. Fig. [4J. 

In spite of these redshift-dependent effects, we still find 
that at all redshifts probed in this paper some fraction 
of red galaxies populate very underdense environments. 
While uncertainties in the environment measurements 
will lead to some galaxies at intermediate densities be- 
ing scattered into the low-density third of the overden- 
sity distribution, upper limits on this effect show that it 
cannot account for all of the red galaxies found in low- 
density environments at any redshift, as we now show. 
We proceed by assuming that the observed log 10 (1 + S3) 
distribution of red galaxies matches the shape of the true 
distribution, an assumption that holds so long as the dis- 
tribution is approximately linear over scales comparable 
to measurement errors (as convolving a linear function 
with a Gaussian leaves it unchanged). In actuality, e.g. 
if the log 10 (I + 63) distribution has some cutoff value at 
low overdensity, the observed distribution will have more 
objects below this cutoff than the true distribution; thus 
the observed distribution provides an upper limit to the 
true contamination rate. 

We calculate this limit from the expectation value of 
the number of red objects that have true overdensities 
in the top two-thirds of galaxies but measured over- 
densities in the lowest one-third, minus the number of 
objects truly in the bottom one-third but measured in 
the top two-thirds. We assume that each overdensity 
measurement may be represented by a Gaussian distri- 
bution in log 10 (1 +^3), with mean given by the mea- 
sured value and standard deviation increasing linearly 
from a = 0.49 dex at z = 0.8 to 0.59 dex at z = 1.2, 
based upon tests with the mock catalogs described in 



§2.4. We find that the increasing uncertainties in envi- 
ronment measures combined with the weakening in the 
strength of the color-density relation with lookback time 
yield a net contamination rate that is roughly constant 
with redshift (~ 35%). Hence, the true red fraction in 
the extreme low-density third of the sample is at least 
3% a t all redshifts. Similar studies at l ower redshift 
(e.g.. [Balogh et al.ll2004bt lYee et a~fll2005l; iG roton et all 
l2005t Martinez et al.ll2006l ) have shown that the existence 
of red galaxies in low-density environments persists to 
z ~ 0. 

Clearly, a population of red galaxies in low-density 
environments exists at 0.75 < z < 1.3. However, 
galaxies may appear red either because they are true 
red-sequence/early-type galaxies or because of the pres- 
ence of interstel lar dust. Based on past studies (e.g., 
lLotz et al.ll2006D . we might expect the former popula- 
tion to dominate the red sequence at lower redshift and 
the latter at higher z. Are the red-sequence members 
in low-density regions at low and high redshift compa- 
rable in terms of galaxy morphology? Out of the less 
than 40 red galaxies in underdense environments identi- 
fied in DEE P2 for which HST / ACS imaging is available 
m the EGS (D avis et al.ll2006f) . an initial by -eye inspec- 
tion indicates that the majority of this population ex- 
hibits early-type morphologies, with the remainder being 
reddened disk galaxies. While the sample imaged with 
HST/ ACS is small, we find no significant trends of mor- 
phology with redshift; for example, dusty disk galaxies 
do not increasingly dominate the red sequence in low- 
density environments at higher z. However, contami- 
nation by late-type galaxies preferentially at higher z 
cannot be excluded as a possible explanation for the ob- 
servational resul ts presented here; an analysis of galaxy 
morphologies by lLotz et alj (|2006l ) concludes that late- 
type galaxies comprise 30% of the red sequence at z > 1, 
though they a re only a minor co ntaminant at low red- 
shift (see also lWeiner et al.ll2005D . Such an increase m 
the contribution of dusty disks to the red sequence at 
z > 1 could work to reduce the strength of an existing 
color-density relation. 

Are members of the red sequence in underdense regions 
the result of passive evolution, in which they consumed 
their gas supplies independent of environmental effects, 
or are they fossil groups that result from the merging 
of several smaller galaxies? Clearly, even in underdense 
environments galaxy mergers are capable of triggering 
events that reduce gas reservoirs and at least temporar- 
ily halt star formation as well as disrupt galactic disks, 
yielding merger remnants with surface-brightness pro- 
files and densit y distributions simi l ar to those of early- 
type galaxies (iToomre fc Toomrei 11972c iNavarro et aT| 
Il987t iBarnes & Hernauistl Il992t bopkins et all I2005TT 
Similarly, passive evolution when teamed with processes 
attributed to secular evolution such as bar instabilities 
could also explain the existence of red, morpholo gically 
early -type galaxies in low-density environments (Zhang 
1996). In future work, we hope to explore the mor- 
phologies of red galaxies in low-density regions within 
the EGS in more detail, with a goal of differentiating be- 
tween these two scenarios; galactic bulges formed by bar 
instabilities tend to differ from those built via mergers 
in that the former exhibit more disk-like properties such 
as flatter profiles and residual bars or spiral structure 
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Fig. 6. — For galaxies in Sample C within high— density (solid line) and low-density (dashed line) environments, we plot the red fraction 
(/r) as a function of redshift for galaxies in sliding bins of Az = 0.1. The high— and low-density samples are selected according to the 
extreme thirds of the local overdcnsity (1 + S3) distribution in the given z bin. The grey shaded regions give the 1-cr range of the red 
fractions in each density regime. The color— density relation evolves significantly over the redshift range 0.75 < z < 1.3, with the red 
fraction in high-density environments dropping continuously as z increases, while in low— density regions it remains relatively constant at 
all redshifts probed. At 2 ~ 1.3, the red fraction shows no dependence on environment within the measurement errors. 
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Fig. 7. — For the mock galaxy samples, we plot the red fraction 
(/r) as a function of redshift in high— density (solid line) and low- 
density (dashed line) environments where the high-density and 
low— density samples are selected as the extreme thirds of the en- 
vironment distribution for the given galaxy sample. The shaded 
grey regions in each plot trace the 1— a range of the red fractions 
in each density regime. Both plots (a) and (b) utilize "true" en- 
vironment measures defined using the entire volume-limited mock 
catalog. However, (a) shows /r(z) for a DEEP2-like "observed" 
sample, while (b) gives the same relation for a sample selected 
according to a fixed (that is, redshift— independent) Mg < —20 
absolute— magnitude limit. Like (a), panel (c) uses a DEEP2-like 
"observed" sample. In contrast, though, panel (c) employs environ- 
ments estimated according to the overdensity (1 + 53) as measured 
in the "observed" sample (i.e., using redshift— space galaxy posi- 
tions and the DEEP2— like galaxy sample). The differences between 
(a) and (c) are consistent with a scenario where noise in environ- 
ment measures causes the color-density relation at high redshift to 
be smeared out. However, such observational selection effects can- 
not alone explain the convergence of that we observe at z ~ 1.3 
in the DEEP2 data (cf. Fig. [HJ. For more details regarding the 
specific sample selection criteria and environment measures, refer 
to the text of §3.1.3. 

(jKormendv fc Fisheiil2005f) . 

4. DISCUSSION 

4.1. Comparison with Related Studies 

From our study of 19,464 galaxies in the DEEP2 
Galaxy Redshift Survey, we conclude that the color- 
density relation observed in the local universe is also seen 
at z > 1, with the fraction of red galaxies increasing with 
local galaxy overdensity at essentially all epochs studied 
(0.4 < z < 1.35) and over all environments probed (from 
voids to large groups). At all redshifts, however, there 
still exists a population of red galaxies in underdense en- 



vironments. In addition, we find that the color-density 
relation evolves with redshift, growing weaker with look- 
back time such that at z > 1.3 the relationship may 
be nonexistent within the range of environments probed 
by the DEEP2 survey (i.e., not including massive clus- 
ters). When viewed in conjunction with the results from 
studies of the galaxy luminosity f unction at < z < 1 
(Bel l et alJl200l iFaber et alJl2005h . our findings provide 
direct evidence that the red sequence is built up preferen- 
tially in overdense environments (i.e., galaxy groups, as 
DEEP2 does not sample rich clusters), thereby produc- 
ing the observed increase in the slope of the red fraction 
versus overdensity relation at later time. 

These results are in agreement with the general picture 
painted by studies of the morphology-density relation 
in galaxy clusters at < z < 1. For ex a mple, build- 
ing up on t he work of iDressler et al.l (|1997t ). ISmith et al.l 
(|2005h and iPostman et al.l (|2005[ ) find that the fraction 
of early-type galaxies increases steadily with density in 
cluster environments out to z ~ 1, with the strength 
of the correlation weaker at z ~ 1 than in loca l sam - 
ples. Corresponding work by Ivan Dokkum et all (|2000D 
also finds significant evolution in the morphology-density 
relation within massive clusters at z < 1, with the early- 
type fraction observed to steadily decline with increasing 
redshift. 

Using data from the Cana da-France-Hawaii Tel escope 
Legacy Survey (CFHTLS), INuiiten et all (|2005l ) simi- 
larly find that over the redshift range < z < 1 
both the red fraction and early-type fraction increase 
in high-density regions with decreasing redshift. While 
the fraction of gal a xies w ith early-type morphologies in 
the INuiiten et al.l (|2005l ) sample is constant with z in 
low-density environments, they find that the red fraction 
(u — g > 1) steadily increases in the "field" from z ~ 0.8 
to z ~ 0, in constrast to our results which show no evolu- 
tion in the red fraction in regions of low galaxy density at 
0.75 < z < 1.3. The differen ces between our results and 
those of lNuiiten et al. (2005) could result from the use of 
photometric redsh ifts to determ i ne env ironments in that 
work; as shown bv lCooper et alj (|2005f ). photometric red- 
shifts cannot cleanly discriminate the environments of 
galaxies, as even the smallest photometric -redshift errors 
achieved are much larger than galaxy correlation lengths 
(< 5/i _1 comoving Mpc versus > 15/i _1 comoving Mpc 
for a z = 0.01). 

Our findings also repr oduce the general trend found 
by iCucciati et ail (|2006l ) based on the VVDS survey; 
they, too, find that the color-density relation weakens 
with redshift over the range < z < 1.5. However, 
within their study the fraction of galaxies on the red se- 
quence shows no significant dependence on overdensity 
at z > 0.9. In contrast, we find a highly significant rela- 
tionship between red fraction and environment at z ~ 1, 
even for highly differing subsamples; only at z > 1.3 are 
our results consistent with density independence, as seen 
in Fig. [51 When binning the DEEP2 data in the same 
redshift ranges as that of ICucciati et al.| (|2006f) . the dif- 
ferences between the DEEP2 and VVDS results are read- 
ily apparent (see Fi gure [8] of this work and Figure 6 of 
ICucciati et al.|[2006h . Given this apparent contradiction 
for z > 0.9 and that the work of ICucciati et all (|2006f) 
employs the most analogous data set to that presented 
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here, a more detailed examination is required. 

First of all, because rest-frame galaxy colors (which 
depend on photometry and coarsely on rcdshift) are 
almost entirely independent of the environment mea- 
surements (which depend upon angular positions and 
high-precision redshifts), we consider it unlikely that the 
highly significant correlation between red fraction and 
environment, which we find at z ~ 1, is false. It persists 
in all four samples considered here, which should differ 
from each other by more than sample B should differ 
from the VVDS volume-limited sample. Since there are 
few galaxies in the gap between the red and blue pop- 
ulations, modest differences in the definition of /r have 
minimal effect, as well, so this is unlikely to explain any 
differences. 

In an effort to conduct a more direct comparison be- 
tween the VVDS and DEEP2 results, we attempt to 
replicate as closely as possible those VVDS subsamples 
for which we can establish a volume-complete DEEP2 
analogue and apply to t hem the same red-fra ction def- 
inition as employed by ICucciati et al.l (|2006( ). While 
in this paper we generally define the red fraction, /r, 
according to a color d ivision in U — B y ersus Mb 
color-magnitude space, ICucciati et alj (|2006[ ) utilize a 
luminosity-independent selection in u* — g' . Using the 
CFHT/Megacam u* and g' filter response, quantum ef- 
ficiency, telescope throughput, and atmospheric extinc- 
tion estim ates 9 , the K-correcti on code (kcorrect version 
v4_l_2) of Blanto n~et alJ (120 03). and our CFHT B R, I 
photometry ( Coil et all l2004bh . we compute the rest- 
frame u* — g' colo r for each galaxy in t he DEEP2 spec- 
troscopic sample. ICucciati et al.l (|2006l ) divide samples 
according to Johnson/ Cousins B-band absolute magni- 
tude in the AB system, taking h = 1; this is identical 
to the Mb used throughout this paper (cf. §2.2). We 
are therefore able to place the DEEP2 galaxies in the 
same color-magnitude space as that of ICucciati et alJ 
(2006) and define "red" according to the VVDS defini- 
tion (u* -g' > 1.1). 

In order to match the VVDS samples as closely as pos- 
sible, we have constructed DEEP2 samples covering iden- 
tical redshift regimes (0.6 < z < 0.9 and 0.9 < z < 1.2) 
and absolute-magnitude limits ( Mb < —20, —20 . 5, —21 ) 
as VVDS subsamples studied by ICucciati et alJ (|2006h . 
Because of the differences between the DEEP2 i?-band 
and VVDS /-band selections, it is possible to replicate 
only some of the VVDS subsamples. For z < 0.9, we 
can construct volume-limited DEEP2 samples down to 
Mb = —20, while for z < 1.2 we are able to construct a 
volume-limited sample matching the VVDS Mb < — 21 
data set, and a very nearly volum e-limited samp l e with 
M B < -20.5. Figures 7 and 8 of ICucciati et all (|2006D 
illustrate the red fraction versus overdensity trends for 
the corresponding samples drawn from VVDS. 

As shown in Figure O the DEEP2 results — for 
volume-limited samples at 0.6 < z < 0.9 and 0.9 < 
z < 1.2 — are consistent with the VVDS results for 
similarly-selected data sets. However, the errors on the 
VVDS trends between red fraction and environment are 
significantly larger (> 2x those for DEEP2). To make 
this figure, we map the data points from Figure 7 of 
ICucciati et all (|2006h onto the DEEP2 results by plot- 

9 http:/ /www. cfht.hawaii.edu/Instruments/Filters 



ting the extreme-overdensity VVDS points at the same 
abcissa values as the extreme-overdensity DEEP2 points. 
This is necessary because citetcucciati06 measures envi- 
ronments over much larger scales than we do in DEEP2, 
and density contrasts on larger scales should be smaller. 
However, locally, at least, the same trends are found us- 
ing environments me asured on smaller and larger scales 
(|Blanton et al.ll2006D . 

No scaling, however, is applied to the VVDS red- 
fraction values, as plotted in Figure [5) The difference in 
normalization of /r between DEEP2 and VVDS could 
be due to K-correction errors in either sample, or, con- 
ceivably, due to a possible tendency of VVDS to fail to 
obtain high-confidence redshifts for red galaxies as of- 
ten as blue, especially at z > 1. From this comparison 
to the ICucciati et al.l (|2006l ) results, we conclude that 
the color-density relation at z > 0.9, as observed by 
DEEP2, is generally consistent with the VVDS measure- 
ments. However, due to their larger uncertainties (likely 
principally due to their smaller sample size), the color- 
density relation found here cou l d not have been detected 
significantly by ICucciati et alJ (|2006t ) . 

Given that a perfect match between our sa mples and 
techniques and those of ICucciati et al.1 ((2006) is impos- 
sible, especially due to the overdensity mapping applied 
in Figure [H we also consider other possible reasons for 
why a trend present in DEEP2 could be missed within 
the VVDS data set. It is quite possible that fundamen- 
tal differences in the samples observed could contribute 
to the differing conclusions at z ~ 1 from the DEEP2 
and VVDS analyses. In particular, the VVDS sample 
includes a subset of objects to Jab = 24, while DEEP2 
samples down to a limit of Rab = 24.1. As a conse- 
quence, while the surveys have comparable depths for 
the bluest objects, VVDS is deeper for red galaxies at 
z > 1, which have (R — I)ab ~ 1-1.5. As a result, 
ICucciati et alJ (|2006h include fainter red galaxies in their 
/r measurements than we do in our standard samples 
(in the volume-limited subsamples, the respective depths 
are indentical, of course), and are able to study samples 
with color-independent Mb limits down to —19.5 (for 
z < 1.2). As shown in Fig. [9l however, volume-limited 
DEEP2 samples matching VVDS samples in luminosity 
range do exhibit a significant color-density trend where 
the equivalent WP S samples do no t. 

In addition to the ICucciati et al.l (|2006| ) sample being 
small (6,582 galaxies) in comparison to our DEEP2 sam- 
ple, yielding increased Poisson errors, it is also limited 
to a single 0.7° x 0.7° field, which increases the uncer- 
tainties due to sample (or "cosmic") variance compared 
to our s ample. The diff e rence s between our results and 
those of ICucciati et al.l (|2006l ), however, are likely not 
attributable to cosmic variance alone. As shown in Fig- 
ure I10[ the variation in the /R-overdensity relation from 
field to field within the DEEP2 survey is small, with each 
field yielding a color-density relation consistent with that 
of the full sample. Cosmic variance mostly changes the 
overall red fraction in a given field and the relative abun- 
dance of the different environments, but not the strength 
of a trend. 

In addition to smaller sample sizes, which yield larger 
errors in /r, there are other phenomena which could 
obscure a true correlation between red fraction and 
overdensity in VVDS data preferentially at high red- 
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Fig. 8. — As in Figure[3] except we plot the red fraction as a function of overdensity and the rest-frame color-magnitude distribution in 
only t wo distinct redshift ranges for t he full galaxy sample (Sample A). Here, the redshift bins are selected to match those of Cucci ati et al.1 
(2006). In comparison to Figure 6 of lCucciati et all (| 2006f ) . we see a highly significant environment trend where they do not. 
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shift. One possibilit y is that environment errors in the 
ICucciati et al.l (|2006f ) sample are large compared to their 
environment bins at z > 0.9. We note that the 5%-95% 
range of measured overdensit ies increases wi t h z in their 
data set, as seen in Fig. 6 of ICucciati et al.l {2006); this 
could reflect increasing errors in their environment mea- 
sures (in the DEEP2 data set, in contrast, the 5%-95% 
range in measured overdensities is nearly independent of 
redshift). 

Another difference is that ICucciati et a l. (2006) use a 
large smoothing kernel (~5 h~ 1 Mpc) in measuring envi- 
ronments. Estimating overdensities on such large scales 
can be problematic due to the small number of inde- 
pendent resolution elements across such a small field (at 
z ~ 1, 0.7° corresponds to ~ 30 h^ 1 comoving Mpc) 
and due to the large percentage of the sample which falls 
within one smoothing scale length of a survey boundary 
or edge (> 50% of the VVDS field by area at z ~ 1). Dif- 
ferences between the true relationship connecting color 
and env ironment on smaller scales and 5 Mpc scales (as 
used by ICucciati et al1l2006f) could also account for dis- 
crepancies between DEEP2 and VVDS results, but at 
z ~ 0, at least, correlations between galaxy properties 
and environment measurements on ~ 1 and ~ 8 Mpc 
scale s appear to differ si gnificantly only in noise proper- 
ties (jBlanton et al.ll2006h . 

Finall y, another effec t that could dilute the measure- 
ments of lCucciati et al.l (|2006l) results from the fact that 
at z > 1 — the redshift regime w here a lack of a color- 
density relation is found in the ICucciati et al.l (|2006l ) 
data — the WDS sample is dominated by lo wer-quality 
VVDS redshifts (flag = 2. lllbert et all 120051 ). A number 
of tests have found that the flag = 2 redshifts have a non- 
negligible error rate (<~ 20%, O. Ilbert, pr ivate communi- 
cation ; C. Wolf, private communication; iLe Fevre et al.l 
l2005al). which wo u ld ca use sources to be included in 
the ICucciati et al.l (|2006h sample at z > 1 incorrectly, 
with both erroneous K-corrections and erroneous over- 
density measurements. This contamination would dilute 
any true correlation between rest-frame color and envi- 
ronment. In this paper, we use only higher-confidence 
DEEP2 redshifts, which have an overall failure rate be- 
low 2% bas ed up on tests with repeated observations 
(jDavis et alJl200a i. 

Wh ile the color-density results of ICucciati et al.l 
(2006) differ when compared to the trends found in the 
full DEEP2 sample, a comparison of similar samples (cf. 
Fig. EH) illustrates that the DEEP2 and VVDS results are 
mostly in agreement where t hey overlap, given th e sig- 
nificantly larger errors for the ICucciati et all (|2006t ) data 
set, which are likely due to the smaller VVDS sample 
size. The various sources of error discussed above (cf. 
§4.1) could also be at play, though the dominant uncer- 
tainty appears to be statistical in nature. Overall, both 
studies independently find that the color-density relation 
grows weaker with increasing redshift. 

Although trends at 0.9 < z < 1.2 are not signif- 
icant for any of the VVDS samples, at lower redshift 
ICucciati et al.l (|2006l ) conclude that the color-density re- 
lation is stronger for brighter galaxies than faint ones. In 
contrast, we do not observe significant variations in the 
color-density trend for samples differing in Mb- That 
is, over the absolute-magnitude ranges probed in Fig. 
[9l the slope of the /R-overdensity relation exhibits no 



statistically significant dependence on luminosity in our 
measurements. The volume-limited subsamples in Fig. 
[SJ however, are overlapping rather than independent 
in luminosity and only probe a limited range of Mb- 
Thus, the lack of significant luminosity dependence to 
the color-density relation within DEEP2 is no t a strong 
statement. As shown by ICooper et al.l (|2006f ). the de- 
pendence of mean environment on luminosity and color 
is effectively separable at z ~ 1, which implies that the 
dependence of the color-density relation on luminosity 
should not be strong over the luminosity range covered 
by DEEP2. 

Within the context of the larger picture of galaxy evo- 
lution, as discussed in §4.2, a luminosity dependence 
to the evolution of the color-density relation may well 
exist. Comparison of our results and those of cluster- 
ing studies at z ~ 1 (e.g.. ICoil et al.l l2004al fbl) to clus- 
tering measuremen ts of bright galaxies at z ~ 2 (e.g., 
iQuadri eT al. 2006) lead t o simil ar conclusions. As ex- 
plained by ICucciati et al.l (|2006h , such a result can be 
easily understood as a n example of cosmic "downsiz- 
ing" (|Cowie et al.l ll996). where the cessation of star for- 
mation occurs first in high-luminosity galaxies or high- 
mass halos (i.e., high-density environments). However, 
the methods used in this paper do not detect this trend 
within the DEEP2 data set. 

4.2. Implications of the Observed Evolution in fn(z) 

A quite striking result from our work is the differ- 
ence between the evolution of the red fraction in under - 
and overdense environments spanning the redshift range 
0.75 < z < 1.3; while /r in underdense environments 
remains roughly constant with z in our sample, we find 
that the red fraction in dense environments decreases 
with z, such that it is roughly equal to the fraction in 
underdense environments at z ~ 1.3. Studying th e prop - 
erties of group galaxies in DEEP2, iGerke et all (2006) 
find similar evolution in the color-density relation from 
0-7 S z S 1-3, with the fraction of blue galaxies in 
groups becoming comparable to that of the field at 
z ~ 1.3. Given ~ 1 Gyr t o move from the blue cloud 
to the red sequence (e.g., Ivan Dokkum fc Franxl l200"ll : 
iBalogh et alj l2004bl ). if the color-density relation first 
arises at z ~ 1.35, then the physical mechanisms respon- 
sible must start to become operative at z ~ 1.7 within a 
concordance cosmology. 

This observed evolution in the color-density relation 
and the existence of an epoch at which environment - 
dependent quenching initiates are both consistent with 
the current theoretical picture in which the conversion 
of blue galaxies into members of the red sequence oc- 
curs in dark matter halos with mass greater than some 
critical value (~ a few x 10 11 M ). As halos pass 
this mass threshold in such hot-flow/cold-flow accre- 
tion models, the infalling cold gas supply to the cen- 
tral galactic disk is virial shocked and shut off such 
that the gal axy will quickly burn its remaining fuel an d 
redden (e.g.. lBirnboim k De kcl 2003: lKeres et al.ll2005h . 
Wit hin the models oflCroton et al.l (l2006h and others (see 
alsolBower et al.ll2006HCattaneo et alj|2006t iKang et all 
I2006T T low-energy AGN activity is included to suppress 
the cooling of shocked gas and the recommencement of 
star formation. 

At yet higher redshifts (z > 2), even in halos 
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Fig. 9. — As in Figure[3] except here we plot the red fraction as a function of overdensity, log 10 (l + S3), in two distinct redshift ranges, 
0.6 < z < 0.9 and 0.9 < z < 1.2. The green circles, red squares, and blue triangles correspond to volume— limited samples with Mb < —21, 
—20.5, and —20, respectively. The shaded (or filled) regions in eac h panel show the 1— a range of the red fraction as a function of overdensity 
for the corresponding VVDS samples, as plotted in Figure 7 of Cucciati et al. (2006); a key to the hatching patterns used is shown at 
upper left. As discussed in the text, the extreme-overdensity VVDS points have been mapped onto the corresponding extreme points from 
DEEP2, while the VVDS red fraction values remain unmodified. Bottom Row: We plot the rest— frame color— magnitude relation, u* — g' 
versus Mb, for all DEEP2 objects in each redsh ift bin. The division b etween the red sequence and the blue cloud is given by the dashed 
line at u* — g' = 1.1, following the definition of Cucciati et al. (2006). The various magnitude-limited samples, Mb < —21, —20.5, and 
—20, are denoted by the green, red, and blue colored symbols, respecti vely. Comparing like galaxy samples, the DEEP2 and VVDS results 
are generally in agreement, given the significantly larger errors for the ICucciati et a l. (2006) data set. 



Evolution of the Color-Density Relation 



17 



0.4 f 

■ • • All Fields 

: Field 1 

: □ □ Field 2 

■ X X Field 3 
0.3 - A-- -AField 4 




0.0 h .... ■ .... ■ .... i .... i .... ■ .... 1 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
l°gio( 1+,5 3) 



Fig. 10. — Using the full galaxy sample (Sample A), we plot 
the red fraction as a function of overdensity, log 1() (1 + £3), at 
0.8 < z < 1 in each of the four DEEP2 fields. The data points 
give the red fraction as a function of the median overdensity com- 
puted in four distinct bins of log 10 (1 + (S3). The horizontal error 
bars run from the twenty— fifth percentile to the seventy-fifth per- 
centile of the overdensity distribution in each bin. The vertical 
error bars give the 1—a uncertainty on the red fraction within each 
overdensity bin, given by Poisson statistics. The lines show linear- 
regression fits to the data points in each field. In the bottom right 
corner of the plot, we give the number of red galaxies within the 
0.8 < z < 1 redshift range in each field. The field-to— field variance 
in the /r— overdensity correlation is small, illustrating the small un- 
certainty due to cosmic variance. Cosmic variance mostly changes 
the overall red fraction in a given field and the relative abundance 
of the different environments, but not the strength of a trend. 

above this critical threshold, quenching does not oc- 
cur since cooling is effectively able t o remove the pres- 
sure support behind the virial shock (Birnboi m fc Dekell 
2003). At l ater times, the evolution of the halo 
mass function (| Jenkins et al.|[200ll iSpringel et al.ll2005l ) 
in combination with the near redshift - independence 
of the critical mass (iKeres et al.l 120051 : ICroton et al.l 
l2006t ICattaneo et alJ l2006h leads to a continual in- 
crease in the number of halos above the threshold mass. 
Given the correlation bet ween halo mass and environ- 
ment in simulations (e.g., [L emson & Kauffnmn^ 119991 : 
iMaulbetsch et alJl2006t IWetzel et alJl2006t ). this thereby 
predicts an evolution in the color-density relation (i.e., 
becoming stronger) at z < 1.3. In this picture, the color- 
density or morphology-density relation within cluster en- 
vironments likely persists out to higher redshifts, as clus- 
ter galaxies reside in the most massive halos, those first 
to reach the critical mass at z ~ 2. As DEEP2 pri- 
marily samples the more common, less massive groups, 
the color-density relation is still weak at somewhat later 
times, z ~ 1.3. 

In addition to predicting evolution of the color-density 
relation at z < 1.3, these hot-flow/cold-flow accretion 
models also provide an additional mechanism by which 
to explain the presence of red galaxies residing in under- 
dense regions (cf. §3.2). While passive evolution or merg- 
ing of several smaller galaxies to form a fossil group are 
also viable physical mechanisms by which to create such 
galaxies, the tail (at > a few x 10 11 M ) of the halo mass 
function in low-density environments naturally leads to 
the existence of at least some true red-sequence galaxies 



in voids. 

The measured evolution in the color-density relation 
as presented here should be directly related to measure- 
ments of galaxy clustering at high redshift. We predict 
that the clustering of ~ L* galaxies at z > 1.3 should 
depend only weakly on color. Studies of galaxy cluster- 
ing by spectral t y pe or by color at h igh redshift (e.g., 
ICoil et al.ll2004allb1 : lMeneux et al.ll2006l ) have shown that 
blue, star-forming galaxies at z ~ 1 are less strongly clus- 
tered than their red counterparts. Such clustering work, 
however, spans a broad redshift range (extending to red- 
shifts less than unity) and the galaxy samples utilized 
are dominated in number by galaxies at z < 1.2, where 
we still find a significant color-density relation. Unfortu- 
nately, the DEEP2 sample includes a modest number of 
galaxies (~ 5000) at z > 1.2, with very few on the red se- 
quence; we therefore lack the statistical power needed to 
compute correlation strengths for subsamples by galaxy 
color at such redshifts. 

At yet higher redshifts (z > 2) there are indica- 
tions that clustering depends on color. Studies of UV- 
selected galaxies and red, near-IR bright galaxies have 
found significant differences in clustering strengths de- 
pending on the sample selection, with the latter be- 
ing more strongly clustered (e.g.. iDaddi et alJ 120031 : 
lAdelberger et al.ll2005bt iFoucaud et al.1 12006). Further- 
more, within near-IR bright samples, the measured cor- 
relation length depends significantly on apparent color 
(e.g., Rab — if Vega or J — K), such that blue near- 
IR bright galaxies cluster like UV-selected Lyman-break 
galaxies (LBGs), while the red near-IR bright galaxies 
exhib it a correlation length larger by a factor of roughly 
two (|Adelberger et al.ll2005at iQuadri et al.ll2006h . 

The red, near-IR bright samples observed at early 
epochs, however, are likely the antecedents of the rarest, 
most massive red galaxies today and are not represen- 
tative of the progenitors of DEEP2 galaxies; their cor- 
relation length at z > 2 already exceeds that of the 
brightest and reddest DEEP2 galaxy samples. Presum- 
ably, if we were to study similarly extreme samples at 
z ~ 1.3 (e.g., galaxies in massive clusters analogous to 
the ca nd idates studied by iRosati et al.l 12004: Mei et all 
200630; iBremer et all l2006f l. we should find a signifi- 
cant color-density relation. Instead, the progenitors of 
DEEP2 galaxies were bluer and/or fainter at z ~ 2 and 
resembled the LBGs more than those red galaxies which 
are most readily observed at high redshift. 

Thus, our observations of the evolution of the red frac- 
tion in low-density and high-density environments do 
not contradict the current set of clustering measurements 
at higher redshift [z ~ 2) or observations of massive 
clusters at z > 1. Like the cluster galaxies, the red, 
near-IR selected, massive galaxies seen at high z, given 
their observed clustering, are likely to reside in very mas- 
sive halos, which should be the first halos to cross the 
threshold quenching mass. In fact, estimates of their 
star-formation rates indicate that they are active ly form- 
ing stars in large quantities (jDaddi et al J 12004 ): hence, 
they do not appear to generally have been quenched at 
z ~ 2, though more recent observations of near-IR se- 
lected galaxies at z ~ 2 indicate relatively low star- 
formation activity in at least a port ion of the massive 
galaxy populatio n at high redshift (|Labbe et al.l 120051 : 
iKriek et al.l[2006h . 
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5. SUMMARY & CONCLUSIONS 

In this paper, we present a detailed study of the evo- 
lution in the color-density relation at 0.4 < z < 1.35. 
Using a sample of galaxies drawn from the DEEP2 
Galaxy Redshift Survey, we estimate the local overden- 
sity about each galaxy according to the projected 3 rd - 
nearest-neighbor surface density. From this, we measure 
the evolution of the red fraction with environment across 
time. Our principal results are as follows: 

• We find that the color-density relation observed lo- 
cally still exists at z > 1 ; the fraction of galaxies on 
the red sequence increases with local galaxy over- 
density to nearly the redshift limits of the DEEP 2 
survey. 

• At all epochs probed (0.4 < z < 1.3), we find there 
exists a population of red, morphologically early- 
type galaxies residing in the the most underdense 
environments. 

• The color-density relation evolves with redshift, 
growing weaker with lookback time such that at 
z > 1.3 there is no detectable dependence of galaxy 
color on local environment in the DEEP2 sample. 

• Our results support a picture in which the red se- 
quence grew preferentially in dense environments 
(i.e., galaxy groups) at z < 1.5. Clearly, the local 
environment plays an important role in "nurtur- 
ing" galaxies, establishing the existence of corre- 
lations such as the morphology-density and color- 
density relation over cosmic time. The strength of 
evolutionary trends suggests that the correlations 
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